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ABSTRACT 

The amount of mass contained in low-mass objects is investigated anew. Instead of 
using a mass-luminosity relation to convert a luminosity function to a mass function, I 
predict the mass-luminosity relation from assumed mass functions and the luminosity 
functions of Jahreiss & Wielen (1997) and Gould, Bahcall & Flynn (1997). Comparison 
of the resulting mass-luminosity relations with data for binary stars constrains the 
permissible mass functions. If the mass function is assumed to be a power law, the 
best fitting slope lies either side of the critical slope, a = —2, below which the mass in 
low-mass objects is divergent, depending on the luminosity function adopted. If these 
power-law mass functions are truncated at 0.001 M@, the contribution to the local 
density from stars lies between 0.016 and 0.039 A4q pc~ 3 , in conformity with the mass 
density that has been inferred dynamically from stars in the Hipparcos Catalogue. 
If the mass function is generalized from a power law to a low-order polynomial in 
log(M), the mass in stars with M < 0.1 Mq is either negligible or strongly divergent, 
depending on the order of the polynomial adopted. 
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1 INTRODUCTION 

The abundance in galaxies of low mass objects is of fun- 
damental importance because we know that near the Sun 
at least half of the Galaxy's mass density is made up of 
stars fainter than Mv = 10, and uncertainty in the expected 
mass-to-light ratio of the Galactic disk is dominated by un- 
certainty in the number of stars with My ^ 16, which are 
extremely hard to detect despite being intrinsically numer- 
ous. 

The traditional way to determine the density of the 
lowest-mass stars is the use of a wide-field proper-motion 
survey to pick up faint but nearby stars, for which pho- 
tometry and perhaps parallaxes can be obtained. More 
recently, two alternative strategies have become available: 
(i) searches for gravitational microlensing events, and (ii) 
narrow-field surveys to identify extremely red stars. 

Microlensing surveys detect stars through their gravi- 
tational fields rather than their radiation, so microlensing 
surveys should provide a powerful probe of the mass con- 
tained in low-mass stars. Unfortunately, there are two large 
problems. First, to obtain even the projected density of de- 
flectors along lines of sight towards the survey stars, one 
must know how the deflectors are distributed along the line 
of sight. Second, to determine the mass function of the de- 
flectors one requires a model of their kinematics. In conse- 
quence of these difficulties, there is no concensus regarding 
the nature of the deflectors that have caused the observed 
microlensing events towards either the Galactic centre or the 
Magellanic Clouds. 

The work of Gould and collaborators with HST (Gould, 
Bahcall & Flynn, 1987; Gould, Flynn & Bahcall, 1988) 
pushes the strategy of counting extremely red stars to its ul- 



timate form, in which the limiting magnitude of the survey 
becomes extremely faint and the field becomes very narrow. 
Consequently, most objects detected are of low luminosity 
and rather distant (~ 2kpc) (as are objects detected by mi- 
crolensing surveys), and it is a non-trivial task, that involves 
the adoption of a large-scale model of the Galaxy, to infer 
the local luminosity function from the data. Despite these 
difficulties, the luminosity function of the Galactic disk is 
now well determined at My < 13 and is usefully constrained 
down to Mv ~ 19 - see Fig. 1. 

This paper is concerned with the problem of convert- 
ing a known luminosity function into a mass function. The 
conventional procedure involves the adoption of some mass- 
luminosity relation. The mass-luminosity relation for cool 
stars is complex and hard to determine either theoretically 
or observationally, while we expect, a priori, that the mass 
function is simple. Therefore, following Kroupa Tout & 
Gilmore (1990) I assume plausible mass functions and use 
measured luminosity functions to infer mass-luminosity re- 
lations. Comparison of these inferred relations with the data 
for binary stars clarifies the amount of mass that may reside 
in sub-stellar objects. Section 2 explains the problem with 
the conventional procedure, Section 3 applies an alternative 
approach to two possible luminosity functions. Section 4 dis- 
cusses the merits of the two luminosity functions, and the 
implications of the mass functions for the local mass density. 
Section 5 sums up. 



2 THE PROBLEM 

Let the number of stars in some volume of space that have 
masses in the range (M,M + dM) be diV = £(M)dM- 
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Then the number of stars that have absolute magnitudes in 
the range (My, Mv + AMv) is clearly given by 



$(My) = 



AN 
AM V 



AM 



AM V 



(1) 



The mass function £ is conventionally determined by esti- 
mating the luminosity function $ by counting stars, and 
then dividing equation (1) through by | | . The problem 
with this procedure is that | goes to zero as one ap- 

proaches the hydrogen burning limit because there M v rises 
extremely rapidly with decreasing mass. Consequently, even 
the smallest error in the assumed value of | will give 

rise to a large error in the derived value of £{M). 

This problem is exacerbated by the fact that Mv{M) 
is not expected to be a smooth function: the atmospheres 
of cool stars are extremely complex with the result that 
their colours vary erratically with temperature and there- 
fore mass. Moreover, low-luminosity stars take up to a Gyr 
to settle onto the main sequence, so Mv really needs to be 
treated as a function of two variables. These considerations 
imply that any empirical determination of | d d ^ j | , for exam- 
ple by fitting a polynomial to empirically determined values 
of Mv(M), will be unreliable. 

Unfortunately, experts in the theoretical modelling of 
low luminosity stars report that the prospects for deter- 
mining | d d j^ | theoretically are no better. For example 
D'Antona & Mazzitelli (1994) write that theoretical T e g 's, 
especially in the red, are intrinsically ill-determined, and no 
sound observational interpretation critically depending on 
the T e ff 's can be presently performed. Clearly, any theoret- 
ical Mv(M) relation does depends critically on theoretical 
T e fi's, and the derivative of this relation will be even more 
uncertain than the relation itself. 



3 AN ALTERNATIVE APPROACH 

The mass function, £(A4) is the outcome of a chaotic pro- 
cess, which involves a wide range of densities, tempera- 
tures, velocities and magnetic field strengths. Consequently, 
it is much more likely to be featureless than either the 
luminosity- mass relation, Mv(M), or the luminosity func- 
tion, $(Mv), which is determined by £{M) and M V {M). 
Hence, it makes sense to assume a simple functional form 
for £(M) and then to use the observed luminosity function 
and equation (1) to predict My{M) and to compare this 
prediction to the relevant observational data. 

From equation (1) we have (Kroupa et al. 1990) 
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(2) 



With this equation I use a generalization of the traditional 
power-law mass function, namely 



= exp fy^Qn^"), 

where 

H = In (M/M e ). 



(3a) 



(3b) 



For N = 1 and en = —2.35 this coincides with the Salpeter 
mass function. For N > 1 the mass function is not a simple 
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Figure 1. The luminosity function of Jahreiss & Wielen (1997) 
(triangles) together with a cubic spline fit (full curve) . Also shown 
are the luminosity functions of Reid et al. (1995), Gould et al. 
(1997), and Gould et al. (1998). The points of Gould et al. (1997) 
have been shifted by 0.1 mag to the left for clarity and joined by 
a cubic spline (dashed curve). 
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Figure 2. The points show empirical determinations of the 
masses of binary components from Popper (1980) and Henry & 
McCarthy (1993). The dashed curve is from Brewer et al. (1993) 
and indicates the predicted locus of subdwarfs. The full curves 
show the mass-luminosity curve one obtains by assuming that 
§ oc M~ 21 . The upper curve is for the Jahreiss & Wielen lu- 
minosity function, while the lower curve is for the Gould et al. 
(1997) luminosity function. The dotted curves are obtained with 
the two luminosity functions and £ oc M~ 21 . 

power law and a characteristic mass is implied by each of 
the coefficients a n for n > 1. 

Fig. 1 shows the two luminosity functions that I have 
used: the full curve is a spline fit to the luminosity function 
of Jahreiss & Wielen (1997; hereafter JW) that employs 
Hipparcos parallaxes. An alternative luminosity function is 
shown by the dashed curve, which deviates from the full 
curve to pass through the points of Gould, Bahcall & Flynn 
(1997; hereafter GBF). 

Fig. 2 shows four curves Mv(Ml) that one obtains for 
the power- law case, N = 1, alongside empirical data points 
from Popper (1980) and Henry & McCarthy (1993). The 
two full curves are for «i = —2.1; the upper curve is for 
the JW luminosity function, while the lower full curve is for 
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the GBF luminosity function. The two dotted curves are 
for Qi = —1.8. The normalizations of <3> and the integration 
constants in equation (2) have been chosen to force all curves 
to coincide at M = 1.6 Mq and 0.09 .M©. The dashed 
curve, which is from Brewer et al. (1993), is the result of 
fitting stellar models to observed subdwarfs. Only one of 
these four curves can be said to be an inadequate fit to the 
data, namely the lower full curve, which is for ai = —2.1 
and the GBF luminosity function. In particular, for My < 4 
the upper full curve (ai = —2.1 and the JW luminosity 
function) fits the data better than the upper dotted curve, 
that assumes the same luminosity function and a = —1.8. 

Formally, none of the curves shown in Fig. 2 provides 
a satisfactory fit to the data - the lowest value of \ 2 per 
degree of freedom is 10.2, which is attained for ai = —2.1 
with the JW luminosity function. Probably no plausibly 
smooth curve would give an acceptable fit to the data, how- 
ever, because at a given value of Mv there are observational 
points that lie far outside the error bars of other observa- 
tional points, probably because there is in reality no single 
mass-luminosity relation: age and metallicity need to be 
taken into consideration. 

Although there is no compelling case to go beyond the 
simplest case N = 1, it is interesting to see what can be 
achieved by increasing iV to 2 and 3. Fig. 3 shows the result 
of minimizing \ 2 with respect to the a„ in equation (3a) 
in the cases N = 2 (full curve) and N = 3 (dotted curve). 
Both fits are for the JW luminosity function and give signif- 
icantly smaller values of \ 2 P er degre of freedom than any 
power law: \ 2 P er degree of freedom is 4.1, 3.1 for N = 2, 3, 
respectively, while the best-fiting power-law has ct\ — —2.06 
and x 2 P er degree of freedom 9.0. (The best power-law with 
the GBF luminosity function has ai — —1.76 and \ 2 P er 
degree of freedom \ 2 — 13.3.) 

Fig. 4 shows the mass functions that underlie the fits of 
Fig. 3. These are remarkably similar for 0.1 < (M/Mq) < 
3, which is the entire range of masses over which the lu- 
minosity function contains information, and are tangent to 
the best-fitting power law at M = 0.4 A4q. The negative 
curvature of these curves clearly implies that the empirical 
mass-luminosity function is better fitted by a mass function 
that at small masses turns down below the best-fitting power 
law. 



4 DISCUSSION 

The luminosity function at My > 13 is controversial. Fig. 1 
illustrates this point by showing in addition to the JW 
points, values derived by Reid, Hawley & Gizis (1995) from 
a kinematically selected sample and by GBF and Gould, 
Flynn & Bahcall (1998) from photometrically selected sam- 
ples of stars observed with HST. The data of Gould, Flynn 
& Bahcall, which are for spheroid stars, only extend to 
My = 13.5 and are compatible with the JW points. The 
points of Reid et al. and GBF go fainter and agree well down 
to My = 14.5. GBF remark that had they been able to de- 
tect more secondaries in binaries, their data would move up 
towards the Reid et al. points, which have been corrected 
for unseen companions. At Mv = 14.5 the GBF and Reid 
et al. data imply that $ is a factor three smaller than the 
value of JW. If the GBF luminosity function were to be pre- 
ferred to that of JW, Fig. 2 would rule out a simple power 
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Figure 3. The same as Fig. 2 except that the curves show the 
result of minimizing \ 2 with respect to the parameters a n in 
equation (3a) for N = 2 (full curve) and N = 3 (dotted curve). 
Both curves are based on the JW luminosity function. 
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Figure 4. The mass functions that correspond to the two curves 
in Fig. 3. 

law with a slope ct\ < —2 that predicts divergent integrated 
mass in low-mass objects. If the JW luminosity function is 
correct, such values of o?i are not favoured over values such 
as ai = —1.8 that place only finite mass in low-mass bodies. 

The luminosity of a stellar population such as the 
Galactic disk is dominated by stars with zero-age masses 
between 0.9 .M© and 3.M©. Therefore, two stellar popu- 
lations that have the same density of stars at Ai = 2A4q 
and above will have very nearly the same luminosity, regard- 
less of their mass functions. Suppose two such populations 
have power-law mass functions for Ai < 2A4q, one of slope 
ai = —2.1 and the other of slope ai = —1.8, and that in 
both cases the power law is truncated at Ai — 0.001 .M©. 
Then the mass, and hence the mass-to- light ratio T, of the 
population with ai = —2.1 will be larger than that of the 
other population by a factor 2.9. Hence the small differences 
between acceptable curves in Fig. 2 are associated with con- 
siderable differences in T. 

From Table 3.13 of Binney & Merrifield (1998; hereafter 
BM) we take that a 2 M@ main-sequence star has My = 1.9, 
and that at this mass {AMy/dhiM) = 3.5. From Table 3.16 
of BM we take that that there are 4.8 x 10~ 4 pc -3 stars per 
cubic parsec with 1.4 < My < 2.4. Hence, with equation (1) 
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a mass function that is normalized to produce the observed 
density of 2 Mq stars is 

C = 8 - 5xl(r4 (2^y '"M^pc- 3 . (4) 

Integrating this down to 0.001 Mq we find that the lo- 
cal mass density in main-sequence stars is 0.039 Mq pc~ 3 
if a = -2.1 and 0.013 .M© pc~ 3 if a = -1.8. Since Ta- 
ble 3.16 of BM gives the V-band luminosity density near 
the Sun as 0.053 L0pc~ 3 , the corresponding mass-to-light 
ratios are Ty = 0.72 and Tv = 0.25. These num- 
bers are usefully compared with the mass density deter- 
mined by Creze et al (1998) from the Hipparcos catalogue: 
0.076±0.015 Mq pc" 3 . Of this latter value ~ 0.04 .M© pc~ 3 
may be contributed by gas, and ~ 0.015 Mq pc -3 may be 
contained in remnants. Hence the main-sequence density 
obtained with a — —2.1 is on the high side, and that ob- 
tained with a = —1.8 is on the low side, but neither number 
can be ruled out, especially given that the mean density of 
local gas is very uncertain. 



5 CONCLUSIONS 

A priori we expect the mass function to be a smooth func- 
tion, while we have every reason to believe that the mass- 
luminosity relation is far from smooth. Indeed, the whole 
concept of a mass-luminosity relation strictly speaking fails 
for an inhomogeneous population such as that of the so- 
lar neighbourhood because luminosity depends on age and 
metallicity in addition to mass. In these circumstances the 
optimum procedure is to assume a simple functional form 
for the mass function and to determine the parameters in 
the fitting function by calculating the mass-luminosity re- 
lation that is obtained by combining it with the observed 
luminosity function. 

This has been done using binary data from Popper 
(1980) and Henry & McCarthy (1993) and two recent lumi- 
nosity functions: that obtained by JW from Hipparcos par- 
allaxes, and that obtained by GBF from HST observations. 
For both luminosity functions, acceptable fits to the data 
can be obtained with a mass function that is a pure power 
law with exponent near the critical value, —2, at which the 
mass implied at small masses diverges. The JW luminosity 
function favours an exponent ai = —2.06 while the GBF 
function favours ai = —1.76. If one continues a power-law 
mass function that is normalized to yield the observed den- 
sity of ~ 2 Mq stars down to M = 0.001 Mq, one obtains a 
local stellar mass density po = 0.039 Mq pc -3 if ai = —2.1, 
and po = 0.013.Mq pc~ 3 if ai = —1.8. The value inferred 
by Creze et al. (1998) from the dynamics of Hipparcos stars 
lies between these figures depending on the rather uncertain 
value of the local density of the ISM. Hence we cannot at the 
present time exclude the possibility that most stellar mass is 
contained in stars less massive than the cutoff for hydrogen 
burning. 

Significantly better fits to the observed mass-luminosity 
relation can be obtained by using a two- or three-parameter 
mass function rather than a power-law. Such functions pro- 
vide a better fit by virtue of their ability to have a slope that 
decreases between Mq and 0.1 Mq. If a two-parameter 
function is employed, negligible mass is contained in stars 



below the hydrogen- burning limit. By contrast, a three- 
parameter function predicts a huge amount of mass in low- 
mass stars. Consequently, nothing can be securely inferred 
about the amount of mass in such stars from currently avail- 
able luminosity functions. The only real constraint is im- 
posed by the dynamics of nearby stars, and this constraint 
is not at all tight. 
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